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We consider core electron photoemission in a localized system, where there is a charge transfer 
excitation. Examples are transition metal and rare earth compounds, chemisorption systems, and 
high Tc compounds. The system is modelled by three electron levels, one core level and two outer 
levels. In the initital state the core level and one outer level is filled (a spinless two-electron problem) . 
This model system is embedded in a solid state environment, and the implications of our model 
system results for solid state photoemission are disussed. When the core hole is created, the more 
localized outer level (d) is pulled below the less localized level (L). The spectrum has a leading peak 
corresponding to a charge transfer between L and d ("shake-down"), and a satellite corresponding to 
no charge transfer. The model has a Coulomb interaction between these levels and the continuum 
states into which the core electron is emitted. The model is simple enough to allow an exact 
numerical solution, and with a separable potential an analytic solution. Analytic results are also 
obtained in lowest order perturbation theory, and in the high energy limit of the semiclassical 
approximation. We calculate the ratio r(aj) between the weights of the satellite and the main peak 
as a function of the photon energy uj. The transition from the adiabatic to the sudden limit is 
found to take place for quite small kinetic energies of the photoelectron. For such small energies, 
the variation of the dipole matrix elements is substantial and described by the energy scale Ed- 
Without the coupling to the photoelectron, the corresponding ratio ro(aj) shows a smooth turn- 
on of the satellite intensity, due to the turn on of the dipole matrix element. The characteristic 
energy scales are Ed and the satellite excitation energy 5E. When the interaction potential with 
the continuum states is introduced a new energy scale Eg = l/{2E?i,) enters, where Rs is a length 
scale of the interaction (scattering) potential. At threshold there is typically a (weak) constructive 
interference between intrinsic and extrinsic contributions, and the ratio r{Lo) /ro{Lo) is larger than its 
limiting value for large to. The interference becomes small or weakly destructive for photoelectron 
energies of the order Es- For larger photoelectron energies r{uj) /roiu) therefore typically has a weak 
undershoot. If this undershoot is neglected, r(ijj) / ro{uj) reaches its limiting value on the energy scale 
Es for the parameter range considered here. In a "shake-up" scenario, where the two outer levels do 
not cross as the core hole is created, we instead find that r{ui) /ro{Lo) is typically reduced for small 
Lu by interference effects, as in the case of plasmon excitation. For the "shake-down" case, however, 
the results are very difi'erent from those for a simple metal, where plasmons dominate the picture. 
In particular, the adiabatic to sudden transition takes place at much lower energies in the case of a 
localized excitation. The reasons for the differences are briefly discussed. 



I. INTRODUCTION 



X-ray photoemission spectroscopy PES is a useful tool 
for studying the electronic structure of solids. The thf 
retical description of PES is however very complicatec 
and almost all work lias been based on the so-called 
sudden approximation.lju The photoemission spectrum 
is then described by the electron spectral function con- 
voluted by a loss function, describing the transport of the 
emitted electron to the surface. The sudden approxima- 
tion becomes exact in the limit when th&^inetic energy 
of the emitted electron becomes infinitefl^Q In this limit 



we can distinguish between intrinsic satellites, appearing 
in the electron spectral function, and extrinsic satellites, 
appearing in the loss function. For lower kinetic energy, 
this distinction is blurred due to interference effects, and 
the satellite weights are expected to be quite different as 
we approach the opposite limit, the adiabatic limit, of 
low kinetic energy.EI It is then interesting to ask at what 
kinetic energy the sudden approximation becomes accu- 
rate. This issue has been studied extensively for the case 
when the emitted electron couples to plasmons, and it 
has been found that the sudden approximation berppies 
valid only for very large (~ keV) kinetic energies.cr □ 
A semi-classical approach has been found to work ex- 
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ceedingly well for the study of plasmon satellites .ffl In 
such a picture, one may take the emitted electron to move 
as a classical particle away from the region where the hole 
was created. The system then sees the potential from 
both the created electron and hole. Initially the electron 
potential cancels the hole potential, but as the electron 
moves away, the hole potential is gradually switched on. 
The switching-on of the hole potential may lead to the 
creation of excitations. If the kinetic energy of the elec- 
tron is sufficiently large, we can consider the hole poten- 
tial as being switched on instantly, and the creation of 
excitations around the hole then reaches a limiting value, 
the sudden limit. In the semi-classical picture, the per- 
turbation is turned on during a time t — Rq/v, where 
V is the photoelectron velocity and i?o is the range of 
the interaction (scattering) potential between the emit- 
ted electron and the excitations. In this picture we also 
need to determine the relevant time scale Tmax so that 
the sudden limit is reached if r ^ Tmax ot v ^ Ro/T^iax- 
Our analysis within the semi-classical framework shows 
that 1/T„iax is related to the energy 5E of the relevant 
excitation of the system and to the strength V of the 
scattering potential. 

We find a different characteristic energy scale Eg = 
l/{2Rl), where Rg is a characteristic length scale of the 
scattering potential. On dimensional grounds one may 
argue that the adiabatic-sudden approximation takes 
place when the kinetic energy of the emitted electron 
is comparable to Eg. This would differ dramatically 
from the semi-classical approach, where the transition 
takes place for energy of the order l/{EsT^ax): i-^- e.g., 
{5EY /Es or V"^ I Es- Alternatively, and again on dimen- 
sional grounds, one may argue that the sudden approx- 
imation becomes valid when the kinetic energy of the 
emitted electron is much larger than tke energy 5E of 
the relevant excitations of the system,Ll in strong con- 
trast to the two criteria above. This latter criterion is 
however not true in general. B 

For many systems with strong correlations, the core 
level spcjctrum can be understood in a charge transfer 
scenario.cl This is illustrated in Fig. 1 for a Cu com- 
pound, e.g., a Cu halide. In the ground-state, Cu has 
essentially the configuration d^ and all the ligand or- 
bitals are filled. In the presence of a Cu core hole, it 
becomes energetically favorable to transfer an electron 
from a ligand to the d-shell, obtaining a d}^ configura- 
tion on the Cu atom with the core hole. Due to the 
hybridization between the d® and d^" configurations, the 
states are actually mixtures of the two configurations, as 
indicated in Fig. |^. In the photoemission process there 
is a nonzero probability that the outer electron will not 
stay on the ligand, but is transferred to the lower energy 
d-like state. This "shake-down" process corresponds to 
the leading peak in the spectrum, while the process where 
the outer electron stays on the ligand corresponds to the 
satellite. This kind-pf model has been appiied-ito rare 
earth compounds,BM chemisorption systems Jlillla transi- 



tion metal compounds and High Tc compounds.^! 
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FIG. 1. Schematic view of the Cu 3s charge transfer pho- 
toemission. Here a is the Cu 3ci level and fe is a ligand (L) 
valence state for the symmetric case. 

Our simple model allows an accurate numerical calcu- 
lation of the photocurrent either by integrating the time- 
dependent Schrodinger equation or by directly inverting a 
resolvent operator (QM). We also derive analytic results 
with a separable potential. These results are compared 
with the semi-classical theory (SC), and with first order 
perturbation theory (PT). In both these cases we have 
analytic results, which is very useful for understanding 
the physics of the problem. 

The impurity model discussed here differs in certain 
important aspects from a real solid. To start with, for a 
solid we never reach the limit of a pure intrinsic spectrum 
since when the cross section for extrinsic scattering goes 
to zero, the range from which the photoelectrons come 
goes to infinity. In our impurity model, on the other 
hand, the extrinsic scattering approaches zero at high 
kinetic energies. Secondly, for a solid, we discuss excita- 
tions in the continuum and not as here discrete energy 
levels. 

For the coupling to plasmons, the adiabatic-sudden 
transition takes place at large kinetiC|-|energies where 
the SC approximation is very accurate. qB The relevant 
length scale is given by the plasmon wave length A ~ 
27r/(7 and the relevant time by the inverse plasmon fre- 
quency Ljq. Large interference effects are then connected 
with a large phase velocity tOq/q, as discussed e.g. by 
Inglesfield£3 Since long wave-length plasmons play an 
important role these large interference effects for small 
q delay the approach to the sudden limit, which only is 
reached at very high kinetic energies (^^ keV). 

For the localized excitations studied here the relevant 
length scale is much shorter. The SC theory then predicts 
that the transition takes place at correspondingly smaller 
kinetic energies. This is indeed what we find from the 
exact solution of our model. This has two consequences. 
Firstly, the SC treatment itself is not valid at such small 
energies, and we have to rely on QM treatments. Actu- 
ally, although the SC treatment correctly predicts a small 
transition energy, we find it predicts qualitatively wrong 
dependencies on the relevant parameters. Secondly, the 
smaller energy scale means that the energy variation of 
the dipole matrix elements becomes very important. The 
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dipole matrix element grows rapidly on an energy scale 
Ed, which can become very important for the adiabatic- 
sudden transition. 

We study the ratio r{uj) between the weights of the 
satellite and the main peak as a function of the photon 
energy lo for the emission from a 3s level. First we con- 
sider the case when the scattering potential between the 
electron and the target is neglected. We find that the 
corresponding ratio rQ{Lu) strongly depends on the ratio 
between the excitation energy 5E and Ed- If ^E/Ed ^ 1, 
ro(w) approaches its limiting value from below, while it 
has an overshoot if 6E/Ed ^1. In both case the limit 
value is reached for photoelectron energies of the order 
of a few times 5E. 

We then study the effects of the scattering potential by 
focusing on r'(a;)/ro(a;). For small energies there is typi- 
cally an overshoot due to constructive interference in the 
"shake-down" case, contrary to the shake-up case where, 
like for plasmons, r(a;)/ro(w) is reduced by interaction 
effects. This happens on the energy scale Eg. If the scat- 
tering potential is very strong, this overshoot may ex- 
tend to several times Eg- Depending on the parameters 
there may be an undershoot for higher energies, which 
can extend up to quite high energies. The undershoot 
is, however, rather small for the parameters we consider 
here, and should therefore not be very important unless 
we want to calculate the spectrum with a high accuracy. 
For Cu compounds and emission from the 3s core level, 
we find that Ed and Eg are comparable, and the relevant 
energy for r(a;)/ro(u;) is then given by £^ ~ iJ^ ~ Ed- 

We present our model in Sec. II and calculate various 
matrix elements in Sec. III. The sudden approximation 
is described in Sec. IV and exact numerical methods 
are given in Sec. V. The perturbational and the semi- 
classical treatments are presented in Sec. VI. In Sec. VII 
we study the condition for the adiabatic-sudden transi- 
tion qualitatively, using simple analytic matrix elements 
and within the framework of the semi-classical theory. 
The results are discussed in Sec. VIII. 



II. MODEL 

We consider a Hamiltonian Tip describing a model with 
a core level c and two valence levels a and 6, 

"Ho = ^aTia + tbUb + tcUc + UaricUa + UbUcnb 

+ i(4c6+4ca). (1) 

The first two terms give the bare energies of the levels a 
and 6, and the last term the hybridization between them. 
The remaining terms involve the occupation number ric 
of the core level c. In photoemission the core level is 
filled in the initial state, and empty in the final, and ric 
only enters as a constant. It is trivial to diagonalize Tio, 
and one obtains two dressed energies Ea{nc) and Eb{nc) 
for the levels a and 6. In a Cu compound, for instance, c 



may represent the Cu 3s core level, a the Cu 3d valence 
level and b a ligand state. This is schematically illus- 
trated in Fig. In our calculations we almost always 
treat the case when Ea{l) > Eb{l), and Eb{0) > Ea{0). 
The meaning of the levels a and b for different types of 
systems with localized excitations is indicated in Table 
p. The full Hamiltonian also has a one-electron part for 
continuum states, 

T = J2 ^k"k, (2) 

k 

with the energies et — k'^/'2,, and wavefunctions '0k ob- 
tained from a one-electron potential corresponding to 
Tic — 0. The rtk are occupation numbers rik = cl^C]^. 
We use atomic units with e = m = h = I, and thus, 
e.g. energies are in hartrees (27.2 eV). The perturbation 
causing photoemission is 

A - 5^(Mk4c, -t- h.c), (3) 

kc 

where Mk is an optical transition matrix element. We 
take the photoelectron interaction as 

^-EK^k?+"<'-^^2]4ck'. (4) 

kk' 

Here V^"^, is a matrix element of the Coulomb potential 
V^^^r) from the charge density /9j/(r) of the orbital v, 

^kk— Jr^ir)V('^Hr)^^^' (r)dr, V^'^Hr) ^ J y^\dr' . 

It is the potential V which determines the transition from 
the adiabatic to the sudden limit, with = we are in 
the sudden limit. 

The total Hamiltonian is given by 

H^Ho + T + V + A. (5) 

This Hamiltonian has two conserved quantities, 

Uc + rik = 1 and Ua + ni, = 1. (6) 

k 

For simplicity we take the core electron and the d elec- 
tron potentials as equal, = V^, and use the relation 
Ua + Ub — 1 to obtain 

^ = "fcE^kk'4ck' (7) 

kk' 

where 

^kk' ^ / ^k(r)l^.c(r)^k' (r)dr, 
Vgcir) EE - y(")(r) = I dr'^^[p,(r') - p,{r')]. 
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The (scattering) potential T4c(r) describes the change in 
the potential acting on the emitted electron when the 
electron in the target hops from level a to b. 
Dropping a constant we can write Hq as 

Ho = ecUc + [sa + Unc)na + etUb 

+ t{cicb + clca). (8) 

where the Coulomb integral U is given by 

U=Ua-Ub 

= I drdr'p,(r)^^[p4r') - pfc(r')]. (9) 



Since the core level is very localized in space this leads 
to 



U = ~VsciO). 



(10) 



For the different types of systems in Table |, a refers to 
a localized level and b refers to a more extended level. 
For instance, for a copper dihalide compound, a refers 
to a Cu 3d orbital and 5 to a combination of orbitals on 
the ligand sites. For simplicity, we approximate the 
ligand orbitals by a spherical shell with the radius Bq 
where Rq is the average Cu-ligand separation. The po- 
tential from the Cu 3d orbital V^dir) can be considered as 
purely Coulombic at r = Rq. The charge from the spher- 
ical shell gives a constant potential inside the radius Rq, 
and we have 



^'"^ ' " \ r>Ro. 



(11) 



Here £ is a constant chosen to make U — —VsciO), which 
may be thought of as being due to screening by the sur- 
rounding. Since V^diO) ^ ^/Ro, e varies only weakly 
with Rq and is approximately given by e ~ V3(j(0)/J7. 



III. MATRIX ELEMENTS 

To estimate the matrix elements Afk and V^^^^' we must 
approximate the photoelectron wavefunctions tp^ (r) . 
These wavefunctions are calculated from the potential 
of a neutral atom, which further is shifted to make the 
potential zero outside a mufhn-tin radius rmt ■ The states 
are then described by spherical Bcsscl functions outside 
rmt , which are matched to a solution of the atomic po- 
tential inside Tmt- For the energy k"^ /2 we obtain the 
partial wave 

{ aiki^ikir) r<rrnt; 
^^^^"^^ 1 \J^k[cosriikjiikr) - sin riikni{kr)] r > r^t, 

(12) 

where ipikir) is the solution of the radial Schrodinger 
equation for the atomic potential inside the mufhn-tin 



radius, aik is a matching coefficient and rjik a phase shift. 
The normalization is given by 



drr^RUr) = 1, 



(13) 



where R is the radius of a large sphere to which the con- 
tinuum states are normalized. The factor (2/i?)^/^fc is 
due to the normalization and the asymptotic behavior of 
ji {x) for large ai_. 

Slater's rulescil are used to generate the orbitals and 
charge densities, from which the potential V^dif) is calcu- 
lated. This gives the scattering potential Vsc{r), which 
is shown in Fig. We consider photoemission from 

a Cu 3s hole. Due to the dipole selection rules, the 
core electron is then emitted into a continuum state of 
p-symmetry. The matrix elements V^.^.' of the scattering 
potential Vsc{r) are shown in Fig. |2p. 



y,,, = / drr^Rk{r)Vsc{r)R,'ir), 



(14) 



where the muffin-tin radius r,„t is taken as the ionic ra- 
dius of Cu, rmt = 2.6 a.u., and we have dropped the I 
index, since we always consider 1 = 1. The dipole matrix 
element Mk is given by 



A4k ^ ak{ek - ec) / drr tpssiry^Jkir)- 



(15) 



We assume that the core level is deep, and jed much 
larger than the energy difference between the ligand and 
copper levels. We can then take the factor Cfc — Cc as a 
constant, which drops out since we always consider rela- 
tive intensities. The result for is shown as the solid 
line in Fig. These dipole and scattering potential ma- 
trix elements are used in the following numerical calcu- 
lations. Extensive calculations of dipole matrix elements 
for many systems were performed by Yeh and Lindau.tj 
To interpret the results, it is useful to also perform an- 
alytical calculations. For this purpose we need models 
of the matrix elements. Below we consider the limits of 
low and high kinetic energies of the emitted electron. In 
the limit of low kinetic energies, we replace the spherical 
Bessel function by its expansion for small arguments 



1 



(2^ + 1)!! 



mix) = ^{21 - ly.l 



1 



,1+1 



(16) 



and the solution ijjik{r) by its zero energy limit ipiQ{r). 
This leads to 



ta.nr]ik 



a-ik 



{krmtr+' 



I + 1 + C (21 - iy.\{2i + iy.\ 

21 + 1 {krmtY 



R L^-f 1 + ^(2^ + 1)!! 



tl^ioirmt) ' 



(17) 
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FIG. 2. (a) The photoelectron scattering potential Vsc{r) 
given by Eq. ( |ll| ) with respect to r for CuCl2 {Ro = 4.71 a.u. 
and e = 1.96). In the inset, we give the atomic configuration 
of Cu-Cl octahedral (nearly octahedral) cluster in CuCb • (b) 
The diagonal and off-diagonal matrix elements of the scat- 
tering potential multiplied by R. In the figure, fci = 1 a.u., 
k2 — 5 a.u., and fcs = 10 a.u. are taken. 

Due to the matching, the coefficient aj^ contains the 
ratio ^ = rmt'>PiQirmt)/'>Pwirmt)- The value of the coeffi- 
cient therefore depends in an interesting way on the wave 
function ipiQ and its derivative, and if Z -I- 1 -I- ^ is close to 
zero aik blows up. Then the matrix elements of the scat- 
tering potential also blow up and we may expect strong 
deviations from the sudden approximation. In such a 
case the dipole matrix element Mik also becomes large, 
i.e., there is a resonance (77/^ = ^) in the photoemission 
cross section. 

From Eqs.(p^) and (|l^) it follows that in the limit of 
a small k the dipole matrix element is proportional to 
Ofc ^ fc'"*"^. The main peak and the satellites in the pho- 
toemission spectrum correspond to different kinetic ener- 
gies and therefore have dipole matrix elements with dif- 
ferent fc-values. This is important at low energies, while 
for large energies the variation in Mik is generally small 
over a range of the energy difference between the main 
peak and the satellite, and for the ratio of the peaks, the 
dipole matrix elements should then not play a role. For 
simplicity, we assume that the dipole matrix elements be- 
come independent of k for R^k 3> 1, where Rd is some 
typical length scale of the system. For our case {I = 1) 



we use the model (note that any constant factor in Mk 
drops out in our final expressions) 



l + (i?rffc)2 " l + Ek/ Ed' 



(18) 



where Ed = l/(2i?^). Fig. |a compares this model with 
the full calculation for a 3s orbital. We obtain Rd — 1-3 
a.u.. For a Is or 2s orbital the length scale is smaller and 
Rd ~ 1/2 a.u.. While we consider I = 1 the behaviour 
for other I is primarily modified for small energies. 

We next consider the matrix elements V^^' . For small 
values of k and k and for I = 1 



^kk' ~ p 



1 2 



V';o(?'mt) 



{kk' 



(e + 2)= 



rRa 


r + 








Jr,„t 






J.2 



Vscir)r^dr. 



(19) 



For small values of k and k it then follows that V^.^,' ~ 
(kk )2. For large values of k and k the matrix elements 
become very small due to destructive interference be- 
tween the two wave-functions unless k k . If A: = fc the 
matrix elements Vkk approach a constant. These features 
are contained in the model 



Vkk' = 



VR, 



[Rlkk 



R [l + {RskY][l + [Rsk'Y][l 



(20) 



where Rg and Rsd are appropriate length scales and V 
has the energy dimension. We recall that V contains in- 
formation about the coefficient aik defined in Eq.(|l3) and 
therefore about the atomic potential and that V may be- 
come particularly large close to a resonance. The dipole 
matrix element Mk and potential matrix V^.^,' in our sim- 
plified model Eqs. (|l^) and (|^ are given in Fig. || as 
compared with the exact results. 

For large values of k and k , the expression Eq. ( |20| ) 
simplifies to 



VRs 



1 



kk 



R l + Rl^{k-k') 



(21) 



An expression of this type can also be derived by assum- 
ing that the wave functions Rik {r) can be approximated 
by spherical Bessel functions in all of space, and by as- 
suming some shape of Vsd'r), e.g., a linear dependence 
on r 
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IV. SUDDEN APPROXIMATION 




-0.2 




FIG. 3. (a) The dipole matrix element Mk as a function 
of k. The exact result (solid line) is obtained from Eq. ( |l^ ) 
and the simplified model by Eq. ( |l8| ) is also shown (dashed 
line). The appropriate parameter is Rd — 1.30 a.u. (b) The 
matrix elements V^.^,/ of the scattering potential are given for 
k — k , k =1 a.u., and k — 5 a.u.. The solid line is from the 
exact calculation for the model of CuCl2 and the dashed line 
is based on the simplified model (Eq.([2o[)). The parameters 
are V — —0.36 a.u., Rs = 1.77 a.u. and Rsd ~ 1-31 a.u.. 



Vsc{r) = K,(0)(1 - — ). 

Uq 

For large values of k and k we then obtain 



^,c(0)i?o l-cos[(fc-A:')ii'o] 



R 



[{k-k')Ro] 



(22) 



(23) 



For this model we relate Rg ~ Rsd — ^o/3 to the range 
Rq of the potential and V — 3K,c(0)/2. Using this iden- 
tification in Eq. (pT|), leads to the correct average value 
of Vkk and to the correct width in k — k of Vj^f.' . The 



simple form (21), however, neglects the effects of the os- 
cillations of the cos-function in Eq. (^3|) for large values 
oi (k — k )Rq, and it therefore gives a worse representa- 
tion of the linear potential (|2^ ) than the form (po| ) gives 
for the more realistic scattering potential (pT|). Later we 
will find that it is a reasonable approximation to put Rd, 
Rsd and Rs equal to the same value R, and introduce the 
corresponding energy E = 1/{2R^). 



We first discuss the photoemission in the sudden limit, 
i.e., we neglect the scattering potential between the emit- 
ted electron and the target {V = 0). The initial state 
|\l/o) is the ground state of Tio with ric — 1 and given by 

|*o) = -sin0|V'c)|V'a)+cos0|Vc)|V'6), (24) 

where 

tan26l = 2t/(e„ + ;7-e6) (25) 
and the corresponding ground state energy is 

Eo^e, + i(e, + U + eb)- ^Vi^a + U - et)^ + At^. 

(26) 

The final states of the target are given by the two eigen- 
statcs of Ti.0 with Uc = 



= cosip\ipa) - smip\ipb), 
I V'2) = sin (filtpa) + COS (p\^pb), 



(27) 



with 



ta.n2ip^2t/{eb~ea) (28) 
and the corresponding energy eigenvalues Ei and E2 are 



El = ^iea + €b)TSE/2, 



with 



SE = y/{ea - ebf + 4^2 



(29) 



(30) 



being the optical excitation energy of the systenx 
The photocurrent J^iuj) {s — 1,2) is given by,El 

Jfe^(L^) = K^f lAlvfo) \^S{uj - ek + Eo- Es), (31) 

where j^E'^'^) is a final state. According to the sudden 
approximation, it can be written as the final target state 
multiphcd by the photoelectron state, jvl'j'^) = \\ps)\''t'k) ■ 
This gives 



(*f |A|*o> 



MkWs 



■nisk, Ws 



— sin (if + 9) , s = 1 
cos (ip + e) , s ^2 

(32) 



J^(ijj) gives the main line (corresponding to the quasi 
particle line in metal) and J^(<^) the satellite line. The 
schematic picture of the initial and final state for this 
system is given in Fig. ^ Summing the kinetic energy 
distribution of the photoelectron, we obtain the absorp- 
tion spectra Js(aj), 
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k,Ws 



k, = y/2{uj + EQ- E,), 



(33) 



where the threshold energies for Ji{lo) and J2(w) are 
given by Ei — Eq and E2 — Eq(= ujth), respectively. 
We can thus also write fci = + 6E — cuth) and 

^2 — \/'2{uj — LUth)- The factor 1/fc conies from the k- 
summation over a (5-function in energy. For convenience 
we introduce the quantity uj = w — w^^, and thus = <D. 

In the sudden approximation the kinetic energy of the 
emitted electron is large, and we can take fci = k2. The 
ratio Too of the satellite to the main peak intensity then 
is 



Too = hm 



- cot^ {ip + e). 



(34) 



Taking into account the energy dependence of the dipole 
matrix element according to model Eq. (18) as well as 
the factor 1/fc, we obtain 



ro(w) = roo 



Cj + 5E 

1 + (il; + 5E)lEa 
1 + Cj/Ed 



(35) 



We now require that the ratio r'o(w) should reach a frac- 
tion 7 (7 ~ 1) of its limiting value ^0(00) for lj — lo^. 
This gives 



if SE<$:Ed; 



6E 



3 1 

2 , "^"'T"' (36) 

j^/^{Ed/SE)*/^ liSE^Ed. 



This criterion refers to the energy where ro(w) reaches a 
fraction 7 in its rising part, and it does not consider that 
there is a large overshoot for SE/Ed ^1. In this case 
we can instead require that ro{u}) is smaller than 7 w 1 
in its descending part. This gives the condition 



1 



6E 



r 



(37) 



In Fig. ^ we show results for ro(w) over a large range of 
values for SE/E. The figure illustrates that the dipole 
matrix element effect alone makes the sudden approxi- 
mation invalid for small kinetic energies. It is interesting 
that for somewhat larger photon energies rg overshoots. 
The reason is that the matrix elements Mk saturate for 
Ck ^ Ed, while the factor 1/fc in Eq.(|3^) favors the satel- 
lite. For 6E/Ed = 1, the result is rather close to the 
sudden limit for uj/5E ^ 1. Finally, for SE/Ed ^ 1, 
there is a substantial overshoot. 
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FIG. 4. The ratio ro(aj) of the satellite to the main 
peak in Eq.(^) divided by the result for an infinite pho- 
ton energy (ro(oo) = roo). Three values of the excitation 
energy SE are considered. 



As discussed in the introduction, we would like to 
study how the adiabatic to sudden transition depends 
on certain factors, like the range Rq of the potential and 
the energy SE of the excitation causing the satellite. We 
therefore keep ratio of t, U and ea—Sb fixed, but vary their 
magnitude. In this way we can vary SE without varying 
the magnitude of the satellite in the sudden limit. Eq. 
( p!o| ) requires that we vary Vsc{0) as we vary SE (via U), 
e.g., by varying the dielectric constant e. In some of the 
calculations below, however, we do not impose Eq. (p^), 
to be able to see the effect of varying SE alone. We fur- 
thermore vary the range Rq of the potential. From the 
definition Eq. (^ij) it follows that this would also vary 
the strength of the potential. For this reason we simul- 
taneously vary the dielectric constant e so that Vsc{0) 
stays unchanged when Rq is changed. Alternatively, we 
can use the analytical matrix elements (^ We can 
then easily vary the length scale by changing R or the 
strength by changing V . 

To know roughly what are interesting values for our 
parameteps, we use experimental results for some copper 
dihalides.E3 We estimate the relative strength of the satel- 
lite to the main peak, and the energy difference between 
the peaks. This gives two equations while in our model 
these quantities, roo and SE, depend on three parame- 
ters, t, U, and Ca — Cb- To only have two parameters we 
consider the symmetric case = e — U /2 and = e as 
shown in Fig.|lj. In the symmetric case we are restricted 
to the shake-down situation since before the transition 
the a-level is above the 6-level, ta + U — eb — U/2, 
while after the transition the a-level is below the h- 
level, €b — €a — U /2. In the symmetric case we have 
< 6* = < 7r/4, and roo = cot^ 2(/3 = U'^/{l&t^). Once 
we know where in the ball-park we have t and U , we can 
leave the symmetric case, and also consider, e.g., shake- 
up cases when there is no level crossing, ea> f-b- In the 
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lowest final state the electron essentially stays on level 
b, while the transfer of the electron to the level a cor- 
responds to a shake up satellite. In this case we have 
-7r/4 < If < < 9 < 7r/4 and ip + 9 <0. 

Our calculations usually take the CuCl2 parameters as 
reference values. For CuCl2 we have 9 — ip = 0.3, which 
gives roo = 2.1. Further V = -0.36 a.u., E = 0.195 
a.u. (with R — 1.6 a.u.), and 6E = 0.237 a.u., i.e., 
V/E = -1.85 and SE/V = -0.66 (see also Table II). 



V. EXACT TREATMENT 

A. Time-dependent formulation 

To obtain exact resu is for model Eq. (^) , we use a time- 
dependent formulationEj and solve the Schrodinger equa- 
tion for the Hamiltonian 

n{T)=no + T + V + Af{T). (38) 

The interaction is switched on at r = 0, using 

/(r)=e-*"^(e-"^-l) ry > 0. (39) 

Here r/ is a small quantity to assure that the external field 
is switched on smoothly. The initial (r = 0) state |^o) 
is given by the ground state of Tio with Uc — 1 in Eq. (^ . 
After a time r, the state |5'(t)) of the system is 

+ ^Cak{T)\i^aMk) +J2^bk{T)\i^bMk)- (40) 
k k 

The coefficients of |^(t)) can be determined by 

,^\^^r))=n{T)\^{r)), (41) 

which gives four differential equations for the four coeffi- 
cients a(r), 6(t), Cafe(T), and Cbk{T), 

d 

^{^a + U + ec)a(r) + i6(r) 

OT 



i-g:^Cbk{T) = (eb + efc)c6fe(T) + tcakir) + V^\T)b{T) 



vY.^f{r)cak{r), 

k 

= {eb + ec)b{T) +ta{T) 
+ Y.^f{r)cbk{r), 



(42) 



(43) 



i-Q^^Cakir) = (ea + efc)cafe(T) + tcbkir) 



Vkk'Cbk'ir), 



(45) 



where Vif{T) = VoMkfir) with Vq representing the 
strength of the external field. We solve the equations 
in the limit when Vq — > 0, and thus the ratio between 
Cak and Cbk is independent of Vq. The initial conditions 
are a(0) = — sin0, 6(0) = cos0, and Cak{0) = Chfc(O) = 0. 
Thus the problem is reduced to solving the coupled dif- 
ferential equations, which is done using the Runge-Kutta 
fourth-order method. 

The photoelectron currents Ji{uj) and J2(w) corre- 
sponding to main and satellite lines, respectively, are 
given by 

Ji(c.)^^|(vl/}'=|*(r))|2 

k 

= I cos (fiCgk (t) - sin (fiCbk (t) | ^ , (46) 

k 

J2(c.)^^|(vl/f|*(r))p 

k 

= ^ I sin (^Cafe(T) cos <^Chfc(r)p, (47) 

k 

where r is a sufficiently large time. We let the system 
evolve for a time of the order 1 /rj to obtain converged re- 
sults for a given finite Vq. In principle, we should perform 
the calculation for a few small values of t] and then ex- 
trapolate to 77 = followed by an extrapolation Vq — > 0. 
Here, for simplicity we have performed the calculation 
for one single small value of Vq. The calculation was per- 
formed for 77 = 0.1 eV, 0.08 eV, 0.02 eV, and the results 
were extrapolated to 77 = assuming the 77 dependence 
a{ijj)rj + b{uj)'r]^ + c{uj). The error in this approach occurs 
primarily for small a;(< 5 eV), and it is then less than 
5% in r{uj)/ro{uj). 

The approach above gives the relative intensity of the 
main and satellite peaks 



r{uj) 



J2M 



(48) 



(44) 



It can be shown that the formulas (^, ^7|) above give 
identical results to the more conventional formulation 
(^ ) below, by performing derivations of the type made 
in, e.g., Ref. |2J. 

As an example of the results obtained in this formal- 
ism, we show in Fig. ^ results for the copper dihalides 
CuBr2, CuCl2 and CuF2. The corresponding param- 
eters are skawn in Table || and were estimated from 
experiment JIj The figure illustrates that there is a small 
"overshoot" for small uj but that the sudden limit is 
reached fairly quickly as Co is further increased. We re- 
mind that in our CuCl2 reference case iJ=0.195 a.u. =5. 3 
eV. 
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(eV) 

FIG. 5. The ratio r(cj) between the satellite and the main 
peak for the divalent copper compounds CuBr2, CuCl2, and 
CuF2. The parameters are given in Table II. The dotted lines 
are the limit values (r(oo)) for the respective cases. 



B. Resolvent formulation 

Alternatively, we can work in the energy space, and 
obtain the spectrum by direct inversion of a resolvent 
operator. We consider the Hamiltonian H, 



H = Ho+T + V, 



(49) 



where by Tip we understand Jig (?ic = 0). The exact final 
photocmission state j^E'j'"') isEl 



1 



E-Ho-T-V-iri 



-V 



mm, (50) 



where (s — 1,2) (Eq. (|27|)) are the exact (target) 
eigenstates of TCq (jif, = 0) and E = Ck + Eg is the energy 
of the final state. Using Eq. (^0|) we calculate the matrix 
element M{s,k) 

M{s, k) 



(*f |A|*o) 



i + v- 



1 



ek + Es-Ho-T -V ^i-n 
Introducing a basis set 

\i) = \i^s)m, 

the matrix elements of V can then be written as 



V, 



Here 



— smiy9, 

C0Siy9, 



if s 
if s 



(51) 
A|*o). 

(52) 
(53) 

(54) 



where we have used Eqs. (|7[ p7|) . The Hamiltonian ma- 
trix in this basis set is diagonalized, which gives the eigen- 
values and the eigenvectors 



We then have 



M(s, k) = M{i) = 
(^|A|*o)+EE 

V 3,1 



(55) 



(56) 



ek + Es - tu + irj 



The quantities (i|A|^'o) were given in Eq. (|3^), 
(i|A|vI/o) = rrii = rUsk = MuWs- By organizing the 
sums in Eq. (|5^) appropriately, the calculation of this 
expression is very fast and the main time is spent in 
diagonalizing the Hamiltonian matrix. We have found 
this method to be more efficient than the time-dependent 
method above. 

In the expression (p6[), we can identify the first term 
as the intrinsic contribution, since this is the amplitude 
which is obtained if there is no interaction between the 
photoelectron and the target. The extrinsic effects are 
then determined by the square of the absolute value of 
the second term. The interference between the intrinsic 
and extrinsic contributions is given by the cross product 
of these terms. 



C. Separable potential 

It is interesting to consider a separable potential 

Vkk'^Vhub^,, (57) 

since it is then possible to obtain an analytical expression 
for r(Lo). The operator in the denominator of Eq. (50) is 
written as 



(58) 



where again \i) — \s)\k) is a combined index for the target 
state s and the continuum state k and z is a (complex) 
number. Then 



and 



dsk{z) = di{z) ^ z - Es 



Csk = bkVs 



(59) 



(60) 



Using the fact that V is separable, it is then straightfor- 
ward to invert the expression in Eq. (|5^) and obtain 



[{z-Ho-T-V)-% 



V- 



(61) 



d^{z) ■ " d,(z)d,(z)(l-UEiC?M(^))" 
This leads to 
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(62) 



1 + sin(p£>fei {Eq + a;)/sin(vp + ^) 
where kg is defined in Eq. 

^,Mfe[l + (l//4)C(e)]' 
with = 1/(2^2) and 

and 



di(e + iri) 



E(e) 



(e + M/) ■ 



(63) 



(64) 



(65) 



To obtain a model for V^.^,' we can, for instance, put 



bk = 




(66) 



Compared with the expression in Eq. (poD, there is no 
termfc — A; in the corresponding expression for V^.^,' . The 
neglect of this term means that Vf.f.i goes to zero more 
slowly as one of the arguments fc or fc goes to infinity. To 
compensate for this we use the power three for Rsk in the 
denominator of (p6|), while in Eq. ( pO| ) the correspond- 
ing power is two. This is a reasonable approximation for 
small k, but it breaks down for large k. 



D. On the variables in the intensity ratio 



For the satellite to main line intensity ratio we have, 



r{uj) 



fci 
k2 



M(2,fe) 



M{l,ki) 



This ratio does not depend on any constant factor in 
Mk, since M(s, k) is proportional to M^- If we take the 
parameters Rd, Rs, and Rsd equal to a common typi- 
cal radius R (as will be motivated later), and use the 
analytic expressions in Eqs. ( |l^ ) and (|o|) then r(w) or 
r(a;)/ro(w), apart from ip and 0, becomes a function of 
SE/E, V/E, and Cj/E, with E = {2lf)-^. We can see 
this since AI^ is a function of kR, and Vkk> a function 
of kR and k'R apart from their prefactors. The prefac- 
tor of Vkk' is VR/R, while that for has no influence. 
For each in a perturbation expansion of Eq. ( |5T| ) we 
have an energy denominator and a fc-summation. The 
/c-summation gives an integral and a factor Rdk. Using 
variables Rk, the R/R in the prefactor vanishes. Factor- 
ing out E in the energy denominator, we have a factor 



V/E for each Vkk', and instead of SE and uj we have 
SE/E and Cj/E. 

With 9 and tf given, SE is proportional to U. U in 
turn is equal to — T4c(0), and thus somehow related to 
the strength of the scattering potential V . If we fix the 
value of the sudden limit roo — cot^ ((^ + 9) by choosing 
one of the angles, we still have an independent parameter 
left. This parameter can be used to decouple the rela- 
tion between SE and V (whatever it is). Summarizing, 
we have found that the parameters of our model system 
appear as the angles 9 and and the excitation energy 
SE (or J7), while the coupling between the photoelectron 
and the model system only appears in one parameter, 
V/E — 2VR^, provided we use Cj/E as variable. We 
have further motivated that we can vary the parameters 
SE and V independently. 



VI. APPROXIMATE TREATMENTS 



A. Perturbation approach to lowest order in V^j 



The same problem can be also studied using the stan- 
dard perturbation approach. We consider the expression 
for the matrix elements M{s,k) in Eq. (^l|). To lowest 
order in V , we can neglect V in the denominator of Eq. 
(pT|). Inserting the completeness relation = 1 in 

terms of eigenstates \i) = \k)\s) we obtain 

M{s,k)^(s\{k\A\^o) (67) 
+ J2^k.,k's' [iE-no-T + iv)-\'s',k's' {s'\{k\A\^o)- 

k a 



Using Eqs. (|53|, |5^) we obtain 

M(l,fc) = ~sm{tp + 9)Mk 
— sin^ (p sm{ip + 9) 



Vkk'Mk' 



sin 2(/3 cos((^ + 8) 



E 

k' 



E - El - Ek' + iri 
Vkk'Mk' 



M (2, k) = cos(.^ -I- 9)Mk 
+ cos^ (fi cos{ip + 9) ^ 

k' 

sin 2(y5 sin((p + 9) ■^-^ 
^ 9 2^ 



E — E2 — Cfc' 



Vkk' Mk' 



IT] 



(68) 



E - E2- Ck' + if], 

Vkk'Mk' 
E - El - €k' + iri 



(69) 



where Vkk' = {k\Vsc\k') and Mj, = {k\A\ipc)- We can 
then immediately calculate the photoemission spectra us- 
ing Eq. 
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.)E= l.VGo\" [f = 0.5 

exact - r[ijj)jrQ[ijj) 

FT - r{ui)!r„[Lo) 

SC - r(u)/roo 



SE = 5.27 c\" {t = 1.3 



SE = 0.14 c\' (i = 1.73 cV) 



iSi3 = 7.02 c\" (f = 2.0 cA") 



<S£ = 7.90 c\' [t = 2.23 c\' ) 



i£ = 8.78 c\" {t = 2.3 



10 20 30 40 
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FIG. 6. r{Lj)/ro{Lo) from the semi-classical approximation 
(SC), the first order perturbation expansion (PT) as well as 
the exact time evolution calculations for different values of 
the excitation energy SE and for U/t = 5.76. The remaining 
parameters are taken from CuCb {Rq^ ~ 4.71 a.u.). 



exact - r(uj)/rQ[uj) - 
PT - riu,)/r„{u,) - 
SC - r(iij)/roo ■ 
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FIG. 7. The same as in Fig. |^ but varying the range DRo 
instead of SE, where D is a scale factor. The parameters of 
CuCl2 are used. 

In Figs. H and |^ we compare the perturbation ex- 
pansion with the exact time dependent calculation for 
a realistic scattering potential in the symmetric case. 



{U/2)y/l + 1/roo, and U = -Vsc{0). Since the ratio 



ro{u}) was discussed extensively in Sec. FV, we here focus 
on r (w)/ro (w) , which describes the effect of the scattering 
potential. We vary the excitation energy SE by varying 
t, while keeping rgo constant. We also vary the poten- 
tial range by replacing i?o in Eq. ( [lT| ) by DRo and then 
varying D. With roo fixed, SE is proportional to Vsc{Q)- 
Thus a small SE and a small D make the perturbation 
weak. The calculations are made for a range of parameter 
values around those given for CuCl2 in Table 



B. Semi-classical approach 

We can also perform the photoemission calcula- 
tion by assuming a classical trajectory of the emitted 
photoelectronEa, producing a time-dependent potential 
which drives the dynamics of the the model. It has been 
reported that the semi-classical approach can give the 
unexpectedlvngcjpd results for the systems with coupling 
to plasmonsJj'llaEa The essence of the semi-classical ap- 
proach is to replace the scattering potential Vsc{r) by a 
time-dependent potential using the charge density p(r, t) 
of the emitted electron, i.e. 



Vscir) 



where we have used p(r, r) 
write the Hamiltonian as 



dvVsc(r)p{r,T)=Vsc{vT), (70) 
S{r ~ vr). We can then 



7i(r) =WoK-0) + y(r), 



(71) 



where Ho^ric — 0) can be expressed in terms of the exact 
final states {tpi and 1P2) in the presence of a core hole 



Hoin, = 0) = Ei^pl^Pi + E2il>lil>2- 
The time-dependent potential takes the form 



(72) 



V{t) - nbVsc{vT) 

= Vll(r)V4V'l + ^^22(t)V'2V'2 + Vl2{T){lp\i:2 + V'lV'l): 

(73) 

where (cf. Eq. (|5|)) 

Vii{t) = Siv? ipVsc{vT), V22{t) = COS^ ipVsc{vT) , 



Vi2{t) = V2i{t) = -]^sin2^Vsc{vT). 



(74) 



The remaining system (target) is still purely quantum 
mechanical, and we write its time-dependent wave func- 
tion |*(t)) as 



|*(t)> ==ai,(T)|Vi)e 



-iEiT 



a2i,(T)|'02)e 



— iE2T 



(75) 



The classical electron velocity v is determined by energy 



conservation, that is, 



1„.2 



We have here chosen 



In the symmetric case we have SE = 2ty^l + roo = the velocity corresponding to the satellite. We could also 
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have performed two calculations, with the velocities cor- 
responding to the leading peak and to the satellite, re- 
spectively. In contrast to the approach used here, this 
would, however, lead to the problem that the spectral 
weight would not be normalized. Applying the time- 
dependent Schrodinger equation to |'I'(t)), we obtain 
aivir) and 021, (r), 

i-^ai^ir) = ViiiT)ai,iT) + l/l2(r)a2„(r)e-'^^^ (76) 



*7r«2.(r) = -^^22^02.^ + y2l(T)al,(T)e^*^^ (77) 

OT 

where 5E (Eg. (|30| )) is the optical excitation energy. 
Eqs. (^6|) and ([ZT]) are subject to the initial conditions 
(cf Eg^) 

0,^(0) = ^,. (78) 

The final photoemission currents Ji{uj) is 

Muj)(x |a,,(To)P, 1 = 1,2. (79) 

It is sufficient to perform the calculation up to t = tq = 
Rq/v, since the potential vanishes for larger values of r. 
The relative intensity between main and satellite contri- 
butions is given by r{uj) — J2{uo) / Ji{uo) as before. 

In Figs. 1^ and |^ we compare the semi-classical and ex- 
act results for a realistic potential in the symmetric case. 
The semi-classical theory is inaccurate over most of the 
energy range considered here. For large energies however, 
the semi-classical theory comes much closer to the exact 
result than does the PT. It is also clear that an increas- 
ing 5E (~ V"sc(0)) does not noticeably affect the energy 
for the adiabatic-sudden transition, where it strongly ef- 
fects the maximum deviation. An increasing D on the 
other hand not only strongly increases the maximum de- 
viation, but also makes the adiabatic-sudden transition 
energy smaller. The dependence on the parameters will 
be investigated more extensively in the next section. 



VII. ADIABATIC-SUDDEN TRANSITION 

We are now in a position to address the adiabatic- 
sudden transition and its dependence on the parameters. 
The calculations are performed with the analytical ma- 
trix elements in Egs. (|l^, |2^). First we study the dif- 
ferent characteristic lengths, 7?^ for the dipole matrix 
elements, and Rs and Rsd for the scattering potential 
matrix elements. We find that it makes sense to use 
only one effective length i?, and the corresponding en- 
ergy E ~ l/(2i?^). As we remarked in Sec. VD, r/r^ 
as a function oi Q/E depends on the parameters 5E/E 
and V /E, and also on the "system" parameters and 
if. We vary 5E independently of V , although for a given 



model there is a direct relation between these two quanti- 
ties. Part of this relation can be offset by using different 
9 and ip (with rpo constant) but we do not explore this 
possibility. The exact solution with a separable potential 
is used to discuss the validity and breakdown of pertur- 
bation theory. We find that V /E has a large effect on the 
deviation from the sudden limit, but little effect on the 
value oi ujjE where the deviation becomes small, while 
5E/E has a comparatively small effect on both magni- 
tude and range of the deviation. For simplicity we use the 
CuCl2 parameters 9 = (p = 0.3, which gives roo = 2.1. 
For CUCI2 we further have V'=-0.36 a.u., £^=0.195 a.u. 
(^=1.6 a.u.), and SE = 0.237 a.u., i.e., V/E = -1.85 
and 5E/V — —0.66. In our calculations, we vary V/E 
and SE/V by typically a factor of two around these ref- 
erence values. 



A. Exact numerical treatment with analytic matrix 
elements 

We first illustrate the dependence on the ratio be- 
tween the length scales Rd, Rs, and Rsd- In Fig. ^ 
we show the results for different ratios Rd/Rs keeping 
Rsd/Rd = 1- These results are obtained for V/Es = —2.0 
and SE/V — —0.5, where Es = l/(2-R?) is the energy 
scale set by the scattering potential length scale. Fig. 
^ shows that as Rd/Rs is reduced the magnitude of the 
"overshoot" is increased. There are, however, no quali- 
tative changes. 



I- 0.8 
0.6 
0.4 





Bd/Bs=1.0 

Bd/Bs=0.8 

Rd/Rs=0-6 




V/E=-2.0 




5E/V=-0.5 



01 2345678 
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FIG. 8. The ratio r{uj)/ro{u;) as a function of Rd/Rs for a 
fixed Rsd/Rs ~ 1 and for (p = 6 — 0.3. The figure illustrates 
that there are no qualitative changes as the length scales for 
the dipole and scattering matrix elements become diflferent. 

Fig. ^ shows results for different values of Rgd/Rs for 
a fixed Rd/Rs — 1. From Eq. ( |20| ) it can been seen that 
this corresponds to varying the range of values k — k 
where Vf.f.i is large, without changing the range over 
which Vkk varies. The figure illustrates that the over- 
shoot becomes larger as Rsd/Rs is reduced. This is nat- 
ural, since decreasing Rsd effectively makes the scatter- 
ing potential stronger by expanding the range of values 
k — k with large scattering matrix elements. The qualita- 
tive behaviour, however, is not changed. In view of Figs. 
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^ and ^, we study below the case when Rsd 
as mentioned in Sec. VD. 



Rs — Rd, 




4 

tt)/E 

FIG. 9. The ratio r{Lj)/ro{uj) as a function of Rsd/Rs for 
a fixed Rd/ Rs = 1 and for = 9 = 0.3. The figure illustrates 
that there is no qualitative changes as the ratio of the two 
length scales in the scattering matrix elements is varied. 

Fig. ^ shows such results for different values of the 
strength of the scattering potential V /E and for differ- 
ent values of the excitation energy 5E/V . In each panel 
5E/V is kept fixed, but the ratio is varied by a factor of 
four from Fig. to Fig. 10:. Typically r(tj)/ro(ijj) has 
an overshoot for small values of Co. For somewhat larger lo 
the ratio approaches unity and possibly becomes smaller 
than unity. The overshoot can be fairly large and hap- 
pens on a small energy scale (~ E). In a few cases of 
a large overshoot, r(aj)/ro(w) does not become approx- 
imately unity until (D is several times E, although the 
relevant energy scale is still E. In the case of an un- 
dershoot, r(coi)/ro(aj) approaches unity from below very 
slowly (energy scale much larger than E). The under- 
shoot is, however, relatively small, and if we do not re- 
quire a high accuracy, we consider the sudden approxi- 
mation is valid when the overshoot becomes small. This 
means that as the range of the scattering potential is 
made larger, the sudden limit is reached at a smaller en- 
ergy. This is the opposite to what one would expect from 
the semi-classical theory. The figure illustrates that 5E 
is not the relevant energy scale. Since in each panel we 
keep SE/V fixed, there is a variation of SE/E by a fac- 
tor of four. Furthermore there is a variation of 6E/V 
by a factor of four in going from the top to the bottom 
panel in Fig. There is no corresponding change in 
the energy for the adiabatic to sudden transition. 



in r(tj)/ro(ixi) quite substantially. Otherwise the results 
are rather similar. For a qualitative discussion, we can 
therefore use the separable potential. 
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FIG. 10. The ratio r{Lj)/ro{ui) as a function of iIi/E for dif- 
ferent values of V/E and 5E/V and ioi ip = 9 = 0.3. The 
figure illustrates that E is an appropriate energy scale for the 
adiabatic to sudden transition. 

For simplicity, we consider 6E = 0. We further put 
Mfc = {Rkf / 1 + {Rkf = hk- This is a poor approxi- 
mation for large fc, but then anyhow also Vkk' is poorly 
represented by the separable potential. Our approxima- 
tions lead to simple results for the functions C, D and E 
entering in Eqs. (|6^|66|). 



k' 



and E{e) = cos9MC{e). Then 



D{e) 



-cosO 



V 



C{e) 



E 1 + {VIE)C{€) 



(80) 



(81) 



Since Dk is independent of k in this approximation, we 
have dropped the index k. The function C(e) is shown 
in Fig. n3. For cp = 6 we can then write 



B. Separable potential 

It is interesting to study a separable potential, since 
it is then possible to obtain an analytical solution. This 
makes it easier to interpret the results. It also allows the 
study the effects of multiple scattering, i.e. the devia- 
tions from first order perturbation theory. Fig. |ll| shows 
results of the exact and first order theory using the same 
values of 6E/V and V/E as in Fig. p^. The separable 
potential overestimates the magnitude of the overshoot 



ro{uj) 



1 



C{u;)] 



1 



2 E 



C{Cu)/[l + ^CiCo)] 



(82) 



For a level ordering as indicated in Fig. |l| < < 7r/4 
and cos^i^/cos(2(/3) > 1. In our standard case with ip = 
0.3, cos^ ip/ cos{2ip) — 1.11. Thus the term in the numer- 
ator of Eq. (H) dominates. The factor 1 -I- {V/E)C{uj) 
gives the multiple scattering, which is not included in the 
first order perturbation theory, i.e. the first order result 
is 
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PT 



cos ^ V 
cos 2ip £; 



1 - 



(83) 



We now compare the behaviors of Eqs. (|82| ) and (|83[ ), 
to see the effects of using perturbation theory. For smaU 
values of w, ReC((D) is positive and then changes sign at 
about Lu/E ^ 2. Im C{lu) is always positive. Due to 
our crude approximations for bk and Mk, C{oj) rapidly 
becomes unreliable beyond oj/E = 2. Both for the exact 
and perturbative expressions r/rQ goes from over- to un- 
dershoot approximately when ReC(tD) — 0. This is some- 
what earlier than in Fig. 11, where however SE = 0.5. 
Comparing r{uj)/ro{uj) for the exact (Eq. (Js^)) and the 
first order result (Eq. (|8^)), we find that the exact solu- 
tion is larger when ReC(a)) > ImC((I'), cf Fig. ^ This is 
consistent with second order perturbation theory, which 
is found to enhance r(uj) /roltu) for small uj and reduce it 
for large uj. For uj = 0, C(Cj) is purely real and slightly 
larger than 0.1, thus multiple scattering gives a diver- 
gence in both Ji(ct') and J2(w) when V/E ~ —10. This 
is due to the attractive potential V forming a bound state 
from the continuum states. 




j/E 



FIG. 11. The ratio r{uj)/ro{u!) for a separable scattering 
potential (^^. According to the exact result in upper panel, 
the overshoot is substantially overestimated by the separable 
potential compared to Fig. px|b, but there is still a qualita- 
tive agreement with the more realistic model (po|). The lower 
panel shows how perturbation theory works and it illustrates 
the effects of multiple scattering. 

The important conclusion from analysing the separable 
potential is that if V is not too large, first order pertur- 
bation theory gives roughly the correct range over which 



there are essential deviations from the sudden limit, while 
multiple scattering increases the magnitude of these de- 
viations for small oj and slightly decreases them for larger 




e/E 



FIG. 12. The function C(e) relevant for a separable po- 
tential is given in the upper panel. The low panel gives the 
behaviors of F{e) defined in Eq.(^). The dotted lines give 
the asymptotic behaviors of F(e) in Eq.(^). 



We next discuss the physical interpretation of the ex- 
pression ( |8^ ) (or the perturbational expressions in ( |67| , 
p^ , IH)). Here unity (the first term in (6^, H], |6^)) cor- 
responds to a direct transition into the final continuum 
state corresponding to energy conservation. The second 
term (the last two terms in ( |67| , |6^ , |69|)) corresponds 
to a virtual transition into some other continuum state 
followed by one or several scattering events with the elec- 
tron ending up in the continuum state corresponding to 
energy conservation. Let us consider the virtual emission 
into a continuum state with a larger energy than the fi- 
nal state and let this be followed by one scattering event 
into the final state. For a negative V the interference with 
the direct event is then constructive. For small photon 
energies such events dominate for two reasons. Firstly, 
there are many more states available above the energy 
corresponding to energy conservation than below, and 
secondly the dipole matrix elements suppress the transi- 
tions to the energies below. As a result, both the main 
peak and the satellite are enhanced by the scattering ef- 
fects. For the values of ip and 9 considered here (< 7r/4), 
the relative effect is stronger for the satellite. As a re- 
sult r(a;)/ro(w) is enhanced. For larger photon energies 
ReC (a)) becomes negative. The density of states of par- 
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tial waves with given / and m quantum numbers decreases 
with energy (~ 1/y/e)). This favours virtual emissions to 
states below the final continuum state. Depending on the 
model for bk and Mk, these matrix elements may have the 
same effect. As a result, Re C((D) becomes slightly nega- 
tive for large energies and the ratio r(Lu) /ro{uj) is slightly 
smaller than one. 

The relevant energy scale for C{uj) is E. This is a com- 
bination of the two effects discussed above. The turn on 
of the dipole matrix elements on an energy scale of the 
order E favors an increasing value of C over this energy 
scale, while the density of states effects becomes more 
important for larger energies. As a result, both Re C{lu) 
and Im C{uj) have a maximum at an energy of the order 
E. 



C. Perturbational treatment with analytic matrix 
elements 



If the more localized level a is above b (see Fig. P in 
the initial state and below b in the final state ("shake- 
down"), we have < (p < ir/A and < 9 < 7r/4, and 
the factor in the brackets is larger than unity. Thus in- 
teraction effects enhance the ratio r{uj)/ro{Lu). This cor- 
responds to a constructive interference between intrinsic 
and extrinsic effects. This is in contrast to-jt-he destructive 
interference found for the plasmon case.DO The present 
treatment, however, refers to the "shake-down" scenario, 
and it is more appropriate to compare the plasmon case 
with the "shake up" case (-7r/4 < ip < < 9 < n/A 
and ip + 9 < 0). Then the expression ( ^7| ) for r{uj)/ro{uj) 
indeed becomes smaller than one, and the relative weight 
of the satellite is reduced for small energies. We notice, 
however, that both the satellite and the main peak are 
enhanced by the interference, but that the main peak is 
enhanced more in the "shake-up" situation. 

We next consider the case when k is large. Fk (e) for 
large k and e is (with Rk « -y/e) 



In this section we study the perturbation theory ex- 
pression in more detail and without relying on a sepa- 
rable potential. Instead we consider the more realistic 
matrix elements in Eqs. (18) and (pO|), assuming that 
Rd = Rs = Rsd = R- We define a function Fk by 



E 

k' 



ZMkFk{e/E), 
E 



(84) 



which is possible due to the simple form of V^^> and M^> 
Explicitly we have 



1 



Rk - 



For the case when 9 — ip we obtain 
r{Lu) \V\1 + 6E/E \V\^ 



ro (t^) 



1 



2a) 2cos(2cp) ACjE 



cos if 1 
cos2(2^ ~ 4 



(89) 



Thus the approach to the sudden limit goes as l/o; with 
a coefhcient which depends on the parameters. With our 
standard CuCl2 parameters, we have 



1 




x'^dx 






1 + 


(^Rk - x^ ^ 


[x^ — e — ii]] 



Fk{e) 



From Eqs. |6^ and |69| we have 
r{uj) 



(85) 



(86) 



1 — cos' if 



2 

E 



(l) 



sin 2cp sin((/)+6') y 



Fk. (^) 



l-sinV|i^fe. {^) 



sin 2ip cos((^4-^) V rp 
2sin(ip+6') "e ^' 



.(§) 



r{uj) \V\ f SE 

-V4 - 1 = -0.30— 1 + 
ro(w) uj \ E 



0.24 



\vl 

CjE 



We note that for large and when \V\ is large 

enough compared to 5E, the last (positive) term dom- 
inates. The approach to the sudden limit is then set by 

\V\''/E = 2[rvY ■ 

To evaluate Eq.(|8^) when 5E = we only need the 
fimction i^(e). 



1 f 

Jo 



l + i^e-xf [. 



Since Rk2 = y Cj/E and Rki = y (a) H- 5E) /E we see, 
as stated earlier, that r/rQ depends only on V/E, SE/E 
and uj/E. 

First we consider the limit of small values of k and 5E. 
For u> = SE = we have F„ (0) = 1/16, and 



We show F (e) in the lower panel of Fig.|ij. The results 
Eq. (|8|) for large values of e are shown by the dotted 
lines in the figure. Clearly the approach of Re F to its 
asymptote is very slow. If we take SE = we have the 
same form as in first order perturbation theory with a 
separable potential Eq. (|83|), 



ro{0) 



1 -I- ^ cost/? cos 0/(16 cos((p -I- 9)) 
1 + ^ sin (yO cos 0/(16 sin(iy9 -I- 9)) 



(87) 



1 - 



COS E 



IV 
2 E 



Fi^) 



(90) 
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where we have put (p = 9. As shown in Fig. F{e) has 
a quaUtatively similar behavior as C(e) for e/E '^2. As 
in the case of the function C, the relevant energy scale is 
E. 
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D. Semi-classical approximation 

In this section we analyze the adiabatic-sudden transi- 
tion within the semi-classical framework. From the cou- 
pled differential equations Eqs. ( fz^ ) and (77) we can ob- 
tain differential equations for d\aii,{T)\'^ /dr, i = 1,2. In- 
tegration of these equations, leads to 



|a2.(ro)|'-|a2.(0)P 



2Im / Vi2{T)aiy{T)a2y{T)e 

'0 



ISEt 



(91) 



where tq ~ Rq/v is the time at which the emitted elec- 
tron with the velocity v leaves the range Rq of the scat- 



tering potential. |aTO(To)|^ 



;(0)p is a measure of the 



deviation from the sudden limit. For small values of r 
both the coefficients aiv{T) and the exponent e^^-^^ are 
approximately real and there is a small contribution to 
the imaginary part of the integral in Eq. (^Tj) . As r grows 
there is, however, a contribution from both these sources. 

To obtain a qualitative understanding of the semi- 
classical approximation, we solve the Schrodinger equa- 
tions Eqs. ( [tg] ) and ( [77|) to lowest order in 1/v. This 
leads to 



FIG. 13. The ratio r{u>)/ro{u!) as a function of lj/E for 
different values of V/E and 5E/V. 

In Fig. |l^ we show r{uj)/ro(uj) as a function of a) /E for 
a few values of SE/V and V/E. We see that r{u})/ro{Lj) 
starts at a positive value (cf Eq. (|7|)), and reaches a 
broad maximum at about uii/E ^ 0.5 — 1.5. Compared 
to the separable potential solution in Fig. the over- 
shoot behavior is robust up to fairly large energies, which 
is due to lmF{uj) decaying more slowly than lmC{uj). 
For much larger values of uj, (r/rp — 1) decays as l/o), 
as shown in Eq. (p9[). Here, based on the discussion in 
Sec. VIIB, we can expect that as multiple scattering be- 
comes important, the region where there is an overshoot 
is substantially reduced and the region with an under- 
shoot becomes larger. At the same time, the overshoot 
intensity will be enhanced. These behaviors are actually 
confirmed by comparing with the exact calculations given 
in Fig. |l^. 

Fig. |l^ illustrates that for intermediate values of u>, 
when RcF dominates, (r/ro — 1) goes as roughly |V^|/a) 
if SE/E is not too large. For larger (but not too large) 
values of ui, ReF becomes small and ImF dominates. 
Since ReF is positive for small energies, this leads to a 
constructive interference between intrinsic and extrinsic 
effects. For energies of the order E, ReF changes sign, 
and the interference becomes weakly destructive. For 
somewhat larger energies the extrinsic effects are ma inly 
determined by the imaginary part of F. From Eq. (p9) 
it follows that in perturbation theory the extrinsic effects 
become small on the energy scale V'^ /E. 



ai.(r) = ai„(0) (92) 

-I f rfT'[l/n(r')ai.(0) + yi2(r')a2.(0)] 
Jo 

and a similar result for a2i,(T). This gives 
|a2.(To)P - |a2.(0)P = ^sin(2(^) |isin(2^) (93) 

/■To "I 2 i-Tq j 

J drVsciT) + sm{2ip + 2e)6E J dTTVsc{T)>. 



To discuss the result, we for a moment assume a simple 
T-dependence of Vij (r) 



To 



(94) 



which corresponds to the r-dependence used in Eq. (|2^). 
We note, however, that this form is too simple to describe 
the behavior of the more realistic potential in Eq. (|lT 
Inserting Eqs. (H) and (||) in Eq. (||) gives 



|a2.(ro)|'-|a2.(0)P 



1 



Vsc{0)sm{2ip) 



i sm{2e)Vsc{0) + i sm{2(p + 2e)5E 



(95) 



We now extend this treatment to intermediate values 
of V where the adiabatic to sudden transition takes place. 
Using the expressions Eqs. (p5|), ( |2^ ) and (30) to relate 
5E and U = — Fsc(O) we obtain 
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r(w) 



- 1 = 



sin(2v?) {SE)'^Tg 



sin2(2v7)cos2(2(^) sin(26l) 
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(96) 



Within the semi-classical theory the condition for the 
sudden approximation is then 



Ro To V 12sin(26') 



(97) 



We are now in a position to discuss the approach to the 
sudden hmit. Within a semi-classical framework it seems 
clear that we have to require that the hole potential is 
fully switched on after a "short" time tq = v/Rq, i.e., 
that the emitted electron leaves the range of the scatter- 
ing potential after a short time. The question is, however, 
what we mean by "short" . From Eq. ( |95| ) it follows that 
the time-scale is set by both the inverse of 6E and the 
inverse of Fsc(O). In Eq. ( |97| ) we have used the relation 
between SE and VsdO) to remove Fsc(O) from Eq. (p7|). 

From Eq. (^) we obtain the condition for the sudden 
approximation within the SC theory 



Q = » {SE)^R, 



{5Ef F,,(0)= 



2E 



E 



(98) 



Thus, according to the semi-classical theory, the sudden 
approximation requires that Cfe ^ {SE)'^/E. Compar- 
ison with the full quantum mechanical calculations in 
Fig. ^ shows that this criterion is not appropriate for 
the range of parameters considered here. The reason is 
that we have considered a parameter range where the 
semi-classical theory is not very accurate. 

It is interesting that the SC theory correctly predicts 
that the weight of the satellite goes to zero at threshold. 
Nevertheless, the SC theory does not give the correct 
physics at the threshold. In the full quantum mechanical 
calculation the weight of the satellite goes to zero due to 
the effects of the dipole matrix element, which becomes 
very small at small photoelectron energies. This effect is 
not included in the SC theory. In the semi-classical treat- 
ment, the small weight of the satellite is due to the fact 
that the scattering potential between the outgoing slow 
electron and the excitation means that the hole potential 
is only switched on slowly. In the quantum mechanical 
treatment, on the other hand, the scattering potential 
leads to an enhancement of the relative weight of the 
satellite close to the threshold for the shake-down case. 



VIII. DISCUSSION 

We have studied the photoemission spectrum of a sim- 
ple model with a localized charge transfer excitation. We 
have obtained exact numerical results for the spectrum 
as a function of the photon energy lu and in particular 
focussed on the ratio r{oj) between the weights of the 



satellite and the main peak. These calculations are com- 
pared with perturbational and semi-classical treatments. 
The results have been analyzed using the latter two ap- 
proaches. 

An important effect in the ratio r{uj) is due to the 
energy dependence of the dipole matrix elements and a 
factor l/{dek/dk) ~ l/fc in the expression for the spec- 
trum. This leads to a suppression of the satellite close to 
the threshold, but can lead to an overshoot further away 
from the threshold. This effect was discussed in Sec. EV 
and is described by ro(w). If the interaction between the 
emitted electron and the target is weak, this effect dom- 
inates. It is determined by the excitation energy SE and 
the relevant energy scale Ed of the dipole matrix element. 
If SE/Ed is small, ro reaches its limiting value from be- 
low, while there is an overshoot if SE/Ed ^ 1. In both 
cases ro reaches its limiting value for a kinetic energy of 
the order a few times SE. 

To study the effects of the scattering potential between 
the emitted electron and the target we have focussed on 
the ratio r(a;)/ro(a;). This quantity shows an overshoot 
for small values of to in the "shake-down" situation stud- 
ied here. Depending on the parameters there may be 
an undershoot for larger energies, which extends over a 
large energy range. This undershoot is, however, fairly 
small for the cases considered here. The sudden approx- 
imation is then valid to a reasonable accuracy when the 
overshoot has become small. We show that this happens 
on the energy scale E ~ l/(2i?^), where i? is a typical 
length scale of the scattering potential. 

One of the main results of this paper is that for a cou- 
pling to localized excitations, the adiabatic to sudden 
transition takes place at quite small kinetic energies of 
the photoelectron. This is in contrast to the large kinetic 
energies needed for the case of coupling to plasmons. In 
the plasmon case, the kinetic energy is typically so large 
that the semi-classical treatment is a very good approx- 
imation. The adiabatic to sudden transition isJiien ex- 
pected to happen on the energy scale {LOqX)'^,M where 
ujq and A are the plasmon frequency and wavelength, re- 
spectively. Since the long wavelength plasmons dominate 
the transition, this happens at very large energies. For a 
localized excitation, the relevant length scale of the scat- 
tering potential is smaller, and the transition is expected 
to take place at a smaller energy scale. Actually, the tran- 
sition takes place at such a small energy that the semi- 
classical theory is usually not valid any more. It is in- 
teresting that the semi-classical theory therefore predicts 
the opposite dependence on the range R of the scattering 
potential, namely as R^ instead oi E = 1/{2R^). 

For the "shake-down" scenario considered here (the 
two outer levels cross as the hole is created), we find con- 
structive interference (increase of r(a;)/ro(w)) between 
the intrinsic and extrinsic processes at low photoelectron 
energies. This is in contrast to the destructive interfer- 
ence found in the plasmon case and to the reduction of 
r(a;)/ro(ijj) found here for the "shake-up" case (no level 
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TABLE I. The model HamiltonianEq.(|l| can describe var- 
ious charge transfer systems. The table indicates the meaning 
of the states a and b for diff'erent cases. 



transition metal compounds 3d state ligand state 

CO on surface 2-k* state bulk(surface) state 

Ce compounds 4/ state 5ci state 

TABLE II. Used parameters for copper-dihalide com- 

tnnds 



pounds 





U (eV) 


t (eV) 


5E (eV) 


Ro (a.u.) 


e 


r{LO co) 


CuBra 


12.33 


2.02 


7.37 


5.01 


1.71 


2.33 


CuCb 


10.58 


1.84 


6.45 


4.71 


1.96 


2.07 


CuFa 


8.62 


1.63 


5.41 


3.86 


2.26 


1.75 
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